IPMU11-0076, NPAC-ll-10,SISSA-48/2011/EP 



Electroweak Beautygenesis: 
From b — > s CP-violation to the Cosmic Baryon Asymmetry 

Tao Liu, 1 Michael J. Ramsey-Musolf, 2 and Jing Shu 3,4 

'Department of Physics, University of California, Santa Barbara, CA 93106, USA 
2 Department of Physics, University of Wisconsin, Madison, WI 53706, USA 
and Kellogg Radiation Laboratory, California Institute of Technology, Pasadena, CA 91125 
''institute for the Physics and Mathematics of the Universe, 
University of Tokyo, Kashiwa, Chiba 277-8568, Japan 
Jt International School for Advanced Studies (SISSA), Via Bonomea 265, 1-34136 Trieste, Italy 

We address the possibility that CP-violation in B s — B s mixing may help explain the origin of 
the cosmic baryon asymmetry. We propose a new baryogenesis mechanism - "Electroweak Beauty- 
genesis" - explicitly showing that these two CP-violating phenomena can be sourced by a common 
CP-phase. As an illustration, we work in the Two-Higgs-Doublet model. Because the relevant CP- 
phase is flavor off-diagonal, this mechanism is less severely constrained by null results of electric 
dipole moment searches than other scenarios. We show how measurements of flavor observables by 
the DO, CDF, and LHCb collaborations test this scenario. 



Introduction The baryon asymmetry of the Universe 
(BAU) has been precisely measured by the WMAP col- 
laboration. Combining its five year results with those 
from other CMB and large scale structure measurements 
gives n b h 2 = 0.02265 ± 0.00059 [I] which is in excel- 
lent agreement with the 95% C.L. range 0.017 - 0.024 
obtained from Big Bang Nucleosynthesis [2 . The im- 
plied ratio of baryon density to entropy s is hb / s ~ 
(8.82 ±0.23) x 1CT 11 . 

To generate the observed BAU, three Sakharov crite- 
ria [3] must be satisfied in the early Universe: (1) baryon 
number violation; (2) C and CP violation; (3) a depar- 
ture from thermal equilibrium (or CPT violation). These 
requirements are not unconquerable, though doing so re- 
quires physics beyond the Standard Model (SM) of parti- 
cle physics. Indeed, there exist a number of possibilities, 
though none have been conclusively established. One of 
the most popular - standard thermal leptogenesis - pro- 
vides a theoretically attractive solution, yet it is generally 
difficult to test experimentally. It is, therefore, worth 
considering scenarios that may be more directly tested 
laboratory experiments. 

A particularly interesting and largely unexplored pos- 
sibility involves CP-violation that enters the B s meson 
system. The relevant phases are generically flavor off- 
diagonal, making them less susceptible to constraints 
from searches for permanent electric dipole moments 
(EDMs) that challenge other baryogenesis scenarios (for 
an illustration in the Minimal Supersymmetric Standard 
Model (MSSM), see, e.g. [I]) . Moreover, recent mea- 
surements in B-factories and at the Tevatron exhibit in- 
dications of CP-violation that differ by a few standard 
deviations from the SM predictions |6], though even 
more recent results from the LHCb collaboration favor 
smaller effects [7J. From a theoretical perspective if the 
CP phase(s) encoding the CP-violation in B s system can 
successfully drive the generation of the BAU and can be 
probed experimentally, our understanding of the BAU 



problem will be considerably advanced. 

In this letter, we report on an initial effort addressing 
this question. We propose a new mechanism in the frame- 
work of Electroweak Baryogenesis (EWBG), explicitly 
showing that the CP-violating phenomena characterized 
by different energy scales (B s observables and BAU) can 
be sourced by a common CP-phase. As an illustration, 
we will work in a Two-Higgs-Doublet Model (2HDM). 
In this context, if a sufficiently strong, first-order elec- 
troweak phase transition (EWPT) occurs in the early 
Universe, the CP phase associated with the tree-level, 
Higgs-6-s interaction at the phase boundary can induce 
CP-asymmetries that ultimately induce the BAU. While 
EWBG in the 2HDM has been discussed extensively in 
the past [H] , including two recent studies using the 2HDM 
that have addressed the possible connection between the 
BAU and B s observables [9] (see also [10]), we empha- 
size that the mechanism discussed below is the only one 
thus far to explore the feasibility of baryogenesis directly 
driven by the b — > s CP- violation. Given the novel fea- 
tures that are generically absent elsewhere and the crucial 
role played by "beauty" quarks, we denote this mecha- 
nism "Electroweak Beautygenesis" (EWBTG). 

In what follows, we concentrate on the issue of CP- 
violation and do not treat the question of the first order 
EWPT in the 2HDM. Following Ref. [9], we instead re- 
fer the reader to more general studies that may indicate 
its feasibility [5j. We note, however, that these anal- 
yses are typically gauge-dependent and therefore open 
to question. Rather than delve into these subtleties of 
perturbative treatments of the EWPT, we also refer the 
readers to a recent and more extensive discussion 
Two Higgs Doublet Model The Higgs sector in the 
general 2HDM is (H u ,d are Higgs doublets with their SM 
gauge charges being (0,2, ±1/2)) 

C = \lQ\eH^ R - \%Q l H d # R 
-yfjQ'HuU^ + yfjQ'ieHDdi + h.c. (1) 
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In a supersymmetric embedding, the first term can arise 
at loop level [12]. For experimental relevance, we focus 
on the two flavor b-s system, with its mass matrix 

rriij = yijV u + XijV d = (y^ sin/3 + A,j cos/3)v, (2) 

where v u 4 are vacuum expectation values (VEVs) of the 
neutral Higgs fields with v — Wv^ + v\ and tan/3 = 
Vu/vd- v u .d are functions of spacetime during the EWPT. 
Meanwhile, Hb s = — cos f3H u + sin/3-ffJ, a linear combi- 
nation of Higgs mass eigenstates, will introduce flavor- 
changing neutral current effects at zero temperature. 

Since we are investigating the feasibility that a com- 
mon phase can source the BAU and account for the B s 
CP-violating observables, we will work in a simplified 
but sufficiently representative scenario, deferring a more 
comprehensive treatment to future work where the fol- 
lowing scenario would arise in one region of parameter 
space. First, we take tan/3 — 1 at zero temperature, em- 
phasizing that tan /3 is not a constant during an EWPT. 
Second, we assume y sb = \ s b = m s b = 0. In the limit 
of y SSl A ss — > 0, there is one CP-violating phase after ap- 
propriate field redefinitions. Without loss of generality, 
we assume that A& s is complex (with 0\ bs = Arg(A{, s )) 
and ybs, ybb and Xbb are real, and furthermore, assume 
\i = Va and \Xb s \ = \Ubs\- The mass matrix is then 



±2£ ss 





±2^6 



(3) 



here £y = \Xij\/V2 and the "±" signs are due to y ss , yb s 
and ybb > or < 0. Denoting rrib s as mb s = Aexp(i#), we 
have A = 2£ bs \cos(6 Xbs /2)\v, = Xb j2 for y bs > 0, and 
A = 2&.|sin(0 Afc ./2)|w, = (0 Xbs + tt)/2 for y bs < 0. 

The mass matrix can be diagonalized by a unitary 
transformation U^MUr = D. In the small m ss limit, 
Ul is trivial and Ur is parametrized by a rotation angle 
cxr = — arctan(A/mf,b). The coupling of Hb s and &l,s,r 
quarks in the mass eigenstate basis is given by 



(bs = 6s[=Fl + exp(z6» Ai ,J] cos Or. 



(4) 



with Arg(Cbs) = #±7r/2 for yb s > and < 0, respectively. 
It is just the phase 8 (or 8 Xbs ) that both introduces CP- 
violation in b — > s transitions (via (bs) and source the 
generation of baryon asymmetry (via rribs)- 
Electroweak Beautygenesis Production of the BAU 
during a first-order EWPT involves a dynamic generation 
of CP-violating charge asymmetries through particle in- 
teractions in the wall of nucleated bubbles. Those charge 
asymmetries are converted, via left-handed fermions 
(n^), into the baryon asymmetry through the elec- 
troweak sphaleron transitions. We ignore the wall cur- 
vature in our analysis so all relevant functions depend 
on the variable z = z + v w t. Here v w is the wall ve- 
locity; z < and > correspond to the unbroken and 
broken phases, respectively; and the boundary extends 



over < z < L w . As pointed out in [15], the trans- 
port properties of particles during the EWPT is most ap- 
propriately treated using non-equilibrium quantum field 
theory. Working in its closed time path formulation (for 
pedagogical discussions, see [5U]) and under the "VEV- 
insertion" approximation (see, e.g., Refs. [51 [T5H2T] ). we 
compute the CP-violating source induced by the Higgs 
mediated process bh — > sr : — > It is given by 

(5) 



U bL LO SR 



xlm- 



(£; L £ SR -e)(n F (£ SR )-n F (£* bL )) 

(£s R -£* br ) 2 



(£b L £ SR + k 2 ){n F (£ SR )+n F {£ bL )) \ 

(£s R +£b L ? J ' 

Here n F (x)=l/(exp(a;) + 1) is the Fermi distribution; 
£b L ,s R — ^b L ,sR — ^b L .s R are complex poles of the spec- 



; and s and 



tral function with w? , -< t »<, . <«n..i .»,,. . . 

DL,S R bL,S R > "L,SR. 

Tb L ,s R are thermal parameters. This source corresponds 
to the "A" -type terms in Eq. (58) of [55J, after properly 
accounting for temperature-independent vacuum contri- 
butions that are removed via normal ordering 31 . The 
quantity 0(z) = d0(z)/dt is given by 



0{z) 



' ( ' Z -sign(y bs )CL(oo) sm0x b 



(6) 



A(z) 2 

with f(z) = (v u (z)v d {z) - v u (z)v d {z)) ~ v w v 2 5(i/L w be- 
ing a function describing the relative variation of the 
Higgs VEVs across the bubble wall. Although analy- 
ses performed in the MSSM [H] indicate 5/3 ~ 0(1CT 2 ), 
a systematic analysis is absent in the 2HDM. Here, we 
will simply adopt 5/3 = —0.05 (if ybs is complex and \bs is 
real, we need (5/3 > to keep the sign of f{z) unchanged). 
Note, Sf^ is non-zero only within the moving bubble wall, 

where 0(z) 7^ 0. 

In contrast to EWBG driven by flavor-diagonal 
sources, the transport of both the second and third fam- 
ily particles is sourced by CP-violating terms. We de- 
fine the number densities {Q\ t 2,3,U,D,C,S,T,B,H = 
£f+ + — — H®} which correspond, respectively, 
to left-chiral quarks of different families, right-chiral up, 
down, charm, strange, top and bottoms, and Higgs 
bosons. Since all light quarks (except &l and sr) are 
mainly produced by strong sphaleron processes and all 
quarks have similar diffusion constants, baryon num- 
ber conservation on time-scales shorter than the inverse 
electroweak sphaleron rate implies the approximate con- 
straints Qi = Q 2 = -2U = -2D = -2C = -2B and 
S + T + Q3 = 0. The set of Boltzmann equations is 



-2r ss (fr - 2£ Q3 + & + 8£b) + Sf L 



T 

m t (£r - (q 3 ) - r t (£r - £ h - (q 3 
S-B) 



-r 



-S^, (with 5 



(7) 



Here 9 M = v w 4= — D a j=2 in the planar bubble wall 
approximation with D a being a diffusion constant, and 
£a = n a /k a with n a and k a being the number density 
and the statistical factor of particle "a" . Apart from the 
CP-violating sources, the interactions in Eq. ^ include 
(i) inelastic top Yukawa (r t ) and strong sphaleron (T ss ) 
processes; (ii) top relaxation processes (r m J, while we 
neglect the other Yukawa interactions since TD q /v^ < 1; 
and (iii) Higgs relaxation processes (T^) due to Higgs 
mass mixing, with typically Th < T mt in this scenario. 

Assuming S^(z < 0) = 0, we solve the Boltz- 
mann equations for the net left-handed fermion density 
iil = Y^i=i Qi analytically order-by-order in 1/T SS , with 
T ss = lQa A s T. The leading contribution arises at first or- 
der in this expansion. The baryon asymmetry p B is then 
produced in weak sphaleron process, described by 



d"PB„ = -e(-2)r„ 
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-pb 



3n; 



(8) 
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FIG. 1: Contours of ub/s (in units of 10~ ). Orange, yel- 
low and green contours indicate 95% C. L. Tevatron (with 
and without di-muon asymmetry) and LHCb constraints, re- 
spectively. We assume kQ 3 — 2kr = 2fcs = 6 |25| . kn = 4, 
v w = 0.4 [2"5]. L.„ = 2/T, D q = 6/T and A bs = 1 TeV. 



where T,, 



120a^T is the weak sphaleron rate 



In 



the broken phase this gives (kg = k B is assumed) 



3T V 



PB 



with r 



D 



dz (9) 



k B (k Q +2k T ) 



and 



AYD- 



_k H (9k T +9kQ+k B ) _ 

v w )/2D. Here D, f and S are, respectively, the effective 
diffusion constant, decay rate and CP-violating source 
for the Higgs number density. D is defined in [5D], while 



f = (9k T + 9k Q + k B )(T mt +2T h )/X 

S = k H (k Q3 -7k T + kB)sffix 
X = 9k Q3 k T + k B kQ 3 + 4k T k B + k B (9kj 



(10) 

9fcg 3 + k B ). 



Note that while the weak sphaleron transitions are driven 
by the diffusion tail for jij that extends ahead of the ad- 
vancing wall in the unbroken phase (z < 0), the solution 
m Eq. (9) contains an integral over the source in the bro- 
ken phase that appears when matching the solutions to 
the Boltzmann equations at the phase boundary. 
CP-violation in the B s — B s Mixing Depending on 
the details of the scalar potential, H bs may be approxi- 
mately a mass eigenstate, which we assume for illustra- 
tion. Tree-level exchange of H^ s with a VEV insertion 
leads to a B s — B s mixing operator in the basis of quark 
mass eigenstates (Af, s is an effective new physics scale): 



Ar 



(b L s R )(b L s R ), with m% 



vA 



bs 



The RG running of this operator involves a mixture of 
scalar operator Og s RR = (bP R s)(bPRs) and tensor op- 
erators 0^ RR = (ba^ v PRs)(ba l 4 V PRs) whose matrix el- 
ements are (B s \O s s RR \B s ) « -hm B J% B h a RR /2A and 



\°trr\ b s) 



-m B Jl Bh s 



I B s TRR 



/2 



and B 



SRR 



r>bs 
TRR 



Assuming m B 
- B B , we ob- 



« (m b + m s 

tain M[ 2 = (B,\H\B.) = -&/§,m B .fli,.(5»?W24+ 
VTRR/2)/A 2 bs , with r) S RR « 1.87, rj TRR « -0.01 [H]. 









0.019 ±0.001 


0.27 ±0.15 


(275 ± 13) MeV 


(Am s J SM 


(Am B j Tev na 


iSM 


19.30 ± 2.2 ps" 1 


17.77 ±0.12 ps^ 1 


(4.2 ±1.4) x 10" 3 


(Ar s ) SM 


(Ar s ) Tcv 


A b 8l m 


0.098 ±0.024 ps" 1 


0.097 ±0.031 ps" 1 


(-7.40 ± 1.93) x 10" 


(Am Ss ) LHCb 


(Ar s ) LHCb 


^LHCb 


17.725 ±0.049 ps" 1 


0.123 ±0.030 ps" 1 


-0.015 ±0.087 



TABLE I: The theoretical input parameters [U and the 
experimental data from the Tevatron and LHCb. 

Choosing T q to be real and parametrizing Mf 2 as [TB] 



(Mf 2 ) SM A s with A s = |A s |e^ , we have 



AT, 



Arf M cos( 
Ar? M sin( 



/SM 



'SM 



Am, = Am 



SM 



A, 



A? 



I A, 



2& = 2/? s bM - </>f (11) 



Here Am s and Ar s are the mass and decay width differ- 
ence between the heavy and light B s mass eigenstates, 
a SL i s the charge asymmetry in semileptonic B s decays, 
and P s measures the time-dependent CP asymmetries in 
the hadronic B s decay. 

The theoretical inputs and experimental results are 
listed in the Table [T] The decay constants and bag pa- 
rameters are taken from Ref. [T5], while A b sl , and (3f xp , 
Ar Exp are obtained by combining the DjZf and CDF mea- 
surements |6]. We perform a % 2 fit to the four ob- 
servables in Eq.(ll|, neglecting the correlation between 
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Ar s and (3 S for simplicity. Assuming Ab s of 1 TeV, we 
scan over the remaining parameters, yielding the regions 
of 95% C. L. from the Tevatron and the LHCb results. 
Discussion The contours of constant tib/s in the 
sign(j/f, s )^h s — sm6\ b3 plane are indicated in Fig. [I] We 
observe that the regions favored by the low-energy fla- 
vor studies at 95% C.L. overlap with regions in which 
a sizable portion of the baryon asymmetry is generated. 
The LHCb results on the B s hadronic decay are more 
constraining on the parameter space than the Tevatron 
ones that do not include the dimuon asymmetry. Al- 
though tension exists between the LHCb- and Tevatron- 
favored regions, it appears feasible that a common CP- 
violating phase may be responsible for both generating 
a non-negligible portion of the BAU and accounting for 
observations in the B s system. 

A definitive statement awaits the resolution of both 
the experimental tensions as well as several theoretical 
issues, including the development of a VEV-resummcd 
CP-violating source (for recent progress, see, e.g., [2"T]). 
analysis of the full numerical solutions to Eqs. Q, and 
completion of a gauge-invariant analysis of the EWPT in 
the 2HDM. Indeed, the results of this initial study are 
likely to indicate the maximum magnitude of the BAU 
that can be achieved in this scenario, given the generous 
assumptions we have made about various input param- 
eters, including Sf3 and v w and the use of an analytic 
rather than numerical solution of the Boltzmann equa- 
tions. Nevertheless, we expect that after future refine- 
ments are implemented, EWBTG may account for an 
interesting portion of the BAU in appropriate regions of 
parameter space. It thus appears that a more compre- 
hensive analysis of this scenario is warranted. 

Though the foregoing discussion relied on the illustra- 
tive case of a two-flavor system of the 2HDM with a sin- 
gle phase, generalization to variants, including minimal 
flavor violation with flavor-blind phases (e.g., see [2"8] ) 
or spontaneous CP-violation, would be straightforward. 
We leave the consideration of these possibilities, along 
with EWBTG in other models such as the four-family 
SM (e.g., see [29]), family non-universal U(l)' model [30] . 
and supersymmetric models, etc. to future work. 
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